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We present a new approach to the problem of Bloch electrons in magnetic field, 
by making explicit a natural relation between magnetic translations and the 
quantum group U q {sl2). The approach allows to express the spectrum and 
the Bloch function as solutions of the Bethe-Ansatz equations typical for com 
pletely integrable quantum systems 



I. INTRODUCTION 

Several times a peculiar problem of Bloch electrons in magnetic field emerged with a new 
face to describe another physical application M, 0, ||, M, []5| 

i) it resembles some properties of the integer Hall effect ||, 

ii) its spectrum has an extremely rich structure of Cantor set, and exibits a multifractal 
behaviour [[|, 0, [0] ( see also || for a review) 

iii) it describes the localization phenomenon in quasiperiodic potential (see e.g. |§ and 
references therein). 

iv) it has been recently conjectured that the symmetry of magnetic group may appear 



dynamically in strongly correlated electronic systems |5), [p~0 1 . 

In this paper we show that this problem (some times called Hofstadter problem) and a 
class of quasiperiodic equations are solvable by the Bethe- Ansatz. We made explicit a long 
time anticipated connection of the group of magnetic translations with the Quantum Group 
U q {sl2) and with Quantum Integrable Systems. The result of the paper is the algebraic Bethe- 
Ansatz equations for the spectrum. Although we do not solve the Bethe-Ansatz equation 



here, we are confident that they provide a basis for analytical study of the multyfractal 
properties of the spectrum. 



II. MODEL AND RESULT 

The Hamiltonian of a particle on a two dimensional square lattice in magnetic field is 

H= £ e^cicR (1) 

<n,m> 

11 e iA *>™ = e l * (2) 

plaquette 

where $ = 2tt^ is a flux per plaquette, P and Q are mutualy prime integers. In the most 
conventional Landau gauge A x = A H - + j = 0, A y = &n x the Bloch wave function is 

^(n)=e^n,(4 ^n=^n + Q (3) 

where n x = n = 1...Q is a coordinate of the magnetic cell. With these substitution the 
Schrodinger equation turns into a famous one-dimensional quasiperiodic difference equation 
("Harper's" equation): 

e ik *ij n+1 + e"^ n _i + 2 cos(k y + n$)^„ = E^ n (4) 

The spectrum of this equation has Q bands and feels the difference between rational and 
irrational numbers - if the flux is irrational, the spectrum is singular continuum - uncountable 
but measure zero set of points (Cantor set). If the flux is rational, then the spectrum has Q 
bands. 

To ease the reference we state the main result of the paper: 

It is known that due to gauge invariance, the energy depends on a single parameter 
A = cos(Qk x ) + cos(Qk y ). We find that the spectrum at A = ( "mid" band spectrum ) is 
given by the sum of roots z\ 

Q-l 

E = tq Q {q-q- 1 )Y J Z h (5) 
l=i 

of the Bethe-Ansatz equations for the Quantum Group U q (sl2) 



i =qQ n tttt.^-g- 1 - ( 6 ) 

with 

q = e^, (7) 

Another version of the Bethe-Ansatz equations is presented in the end of the paper. Solution 
of the model with anisotropic hopping ( when the coefficient in front of cos in Eq.(4) differs 
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from 2 is also available. It will be published elsewhere. The Quantum Group symmetry is 
more transparent in another gauge A x = — ®(n x +n y ), A y = ®(n x + n y + 1). In this gauge 
a discrete coordinate of the Bloch function if)(n) = exp (ipn)ip n (p) , turns to n = n x + n y . 

It is defined in two magnetic cells : n = 1...2Q, p± = (p x ±p y )/2 e [0, $/2]. An equivalent 
form of the Harper Eq. (^) for ip n is 

2e i*+*P+ cOs(i$n + ^$ - P-)iJ n +l + 

2e~i*- ip + cos(i$n - ^ - p_tyn-i = (8) 

In the new gauge the "midband " A = corresponds to the point p = (~7r, |7t) 0. The 
advantage of this gauge is that the wave function turns into the polynomial with roots z m 

^(z)=f[(z-z m ) (9) 

m=l 

at the points z = q 21 

Vv = *(q n ) (io) 



III. GROUP OF MAGNETIC TRANSLATIONS 

The wave function of a particle in a magnetic field forms a represantation of the group 
of magnetic translations |l| : let generators of translations be 

T p (i) = e iA ^+p \ ixi + fl\ (11) 

They form the algebra 

-L ft —fit -'n-'m y -^n+m> 

TyT x = q T x T y , T y T_ x = q T^. x T y ... (12) 

with q given by the Eq.(7). The Hamiltonian ([I]) can be expressed 

H = T X + T_ x + T y + T_ y (13) 

This group is also equivalent to the Heisenberg-Weyl group: [p,q] = T x = expg, T y = 
exp p. 



1 The doubling of period in comparison with original Harper's equation is artificial: using a simple 
transformation (multiplying ip n by e~ l ® n / 4 ) one comes to an equation with coefficients of period 
Q. We are indebted to Alexandre Abanov for clarifying connection between two gauges 
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IV. QUANTUM GROUP 

The algebra U q (sl2) ( a q-deformation of the universalenveloping of the si 2) is generated 
by the elements A,B,C,D, with the commutation relations |lTJ, |TJ], [|13|], |^4| , |T5|| 



AB = qBA, BD = qDB, 
DC = qCD, CA = qAC, 

a 2 - n 2 

AD = l, [3,0] = —^ 

(14) 

The center of this algebra is a q— analog of the Casimir operator 

In the classical limit q — ► 1 + |$, the quantum group turns to the SI2 algebra: 
(A - D)/{q - q- 1 ) -^S 3 ,B^ S+, C^S.,c^S 2 + 1/4. 

The commutation relations ([14]) are simply another way to write the Yang-Baxter equa- 
tion 

where generators A, B, C, D are matrix elements of the L- operator 

- uA-u~ 1 D u~ 1 C ' 

Here u is a spectral parameter and R- matrix is the L- operator in the spin 1/2- 
representation. It is given by the same matrix ( |P7D with elements: A = q2 a:i , D = 
q~2 (73 ) B = cr + , C = cr_, where a are the Pauli matrices. 

Finite dimensional representations (except some representations of dimension Q ) of 
the U q {sl2) can be expressed in the weight basis, where A and D are diagonal matrices: 
A = diag (q\ q~i). An integer or halfmteger j is the spin of the representation, and 2j + 1 
is its dimension. The value of the Casimir operator ( |T5D in this representation is given by 
the q- analog of (j + 1/2) 2 

°=( 3 ^ ) =b- + i/2i; (is) 

Representations can be realized by polynomials ^(z) of the degree 2j: 
A9(z) = q~ j ^(qz), D^(z) = q j '^(q' 1 z), 
BV(z) = z(q - q- 1 )- 1 (q 2j V(q~ l z) - q~ 2j ^(qz)) 

CV(z) = ~z~ 1 {q - q- 1 )- 1 (^(q^z) - V(qzj) (19) 



This again the g-analog of the representation of the 5/2 algebra by a differential operator: 

Ss = z Tz~ jj s+ = z ^ 3 ~ z J^' S - = Tz (20) 
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V. MAGNETIC TRANSLATIONS AS A SPECIAL REPRESENTATION OF THE 

QUANTUM GROUP 



Dimension of our physical space of states is 2j + 1 — Q. This is a very special dimension 
when g 2j+1 = =pl for P - odd (even). The Casimir operator fll8|) in this case is 

c = -4(g - q' 1 )' 2 , for P - odd 

c = 0, for P — even (21) 

In this special case |16| , fl7j representation of the quantum group can be naturally (but not 
unambiguously ||18|| ) expressed in terms of Magnetic Translations. Say one may choose 

T_ x + T_ y = ±i(q - q~ l )B, 
T x + T y = i(q-q~ l )C, 

P-yPx = ±q~ l A\ T- x T y = ±qD 2 (22) 

where the upper sign corresponds to an odd P and the lower to an even P (representation 
of U q (sl2) in terms of different but related Weyl basis can be found in Ref. fl6|| , IPTfl , |I8| , 
f2lfl). It is straightforward to check that this representation obeys commutation relations 
(|l2l) and (|l4l) and gives correct values (21) of the Casimir operator 



The Hamiltonian (f[,|13) now can be expressed in terms of the quantum group generators 

H = i(q - q- l )(C ± B) (23) 

whereas the Schrodinger equation becomes a difference functional equation 

i{z~ x + qz)^(qz) - i{z" x + q~ l z)^(q~ l z) = E^(z) (24) 

The original Harper's discrete equation in the form (|| ) at the point p = (7r/2,7r/2) can 
be obtained from the functional equation (|24| ), by setting z = q l and ipi = ^/(q l ). The 
advantage to use the extention of ipi to a complex plane z is that the representation theory 
of the quantum group garantees that in a proper gauge the extended wave function would 
be polynomial (9). 



In addition to representation (19) having the the highest and the lowest weight , in 



the special dimension g 2j+1 = ±1 there is a parametric family of representations having 



in general no highest or lowest weights | 17| , ||21|| . The parameter describes the anisotropy 



of the hopping amplitude in the Hamiltonian (1) or the strength of the potential (i.e. the 
coefficient in front of cos) in the Harper Equation (4). In this case the wave function is 
not polynomial in any gauge. Nevertheless the Bethe Ansatz solution of the anisotropic 
problem is also possible but acqueres much heavier mathematical techique. We postpone it 
for a more extended paper. 

VI. FUNCTIONAL BETHE- ANSATZ 

Among various methods of the theory of quantum integrable systems the functional 



Bethe-Ansatz [|n| seems to be the most direct way to diagonalize the Hamiltonian (|2 
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We know that the solution of Eq.(|24]) is polynomial 



m{z)=X[{z-z m ) (25) 



m=l 



at the points z = q 21 

V>» = *(g n ) (26) 
Let us substitute it in the Eg . (|24l) and divide both sides by ^(z). We obtain 



Z Z m 



-i(z- l + q- l z) \[ ^ — = E (27) 

m=l,mj^l m 

The l.h.s. of this equation is a meromorphic function, whereas the r.h.s. is a constant. To 
make them equal we must null all residues of the l.h.s.. They appear at z = 0, at z = oo 
and at z = z m . The residue at z = vanishes automatically. 

The residue at z = oo is — iq® + iq~® . Its null determines the deqree of the polynom. 

Comparing the coefficients of z®~ x in the both sides of Eq.(24),we obtainthe energy given 
byEq.(5). 

Finally, annihilation of poles at z = z m gives the Bethe-Ansatz equations (|6]) for roots of 
the polynomial (9). We write them here in a more conventional form by setting z\ = exp (2<pi) 

cosh(2<^-if) 9z} sinh(^ - ip m + if ) 

cosh(2<^ + z T ) w= i, m ^ smh{<pi - (p m - ij) 

Another form of the Bethe-Ansatz equations is given in the next section. 



VII. MISCELLANEOUS RESULTS 
1. Another Form of the Bethe-Ansatz Equations 

As we already mentioned the representation of the quantum group by magnetic trans- 
lations is not unique. This means that there is another gauge where the wave function 
is a polynom. Consider the gauge A x = —A y = -$n x . Then instead of Harper's 
equation^) we obtain an equivalent equation for the gauge transformed wave function 
ipn — > (fin = exp(i^n(n — l))^ (for an odd Q it respects the periodic conditions (3)). 
At A = 0, it has the form 

e-^ n (fi n+1 + e^^Vn-i - 2 cos(ri$)0 n = Ecfi n (29) 

This choice of the gauge corresponds to another representation of the quantum group by 
magnetic translations. Say for an odd P we have 
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T_ x + T_ y = -t{q-q- 1 )q-^BD, 
T x + T v = -i(q-q- 1 )q-$CA, 

T- y T x = q-'A 2 , T_ x T y = qD 2 (30) 
Then the Hamiltonian (jl3| ) turns into quadratic form in the U q (sl2) generators 

H=-i(q-q- 1 )q-*(CA + BD) (31) 

The representation (|T9D (for g 2j+1 = — 1) now gives another functional equation 

z' 1 ^{q 2 z) + q~ 2 z^{cf 2 z) - (z + Z - 1 )m{z) = E^{z) (32) 

which is identical to ( p9|) on the set of points z = q 21 , where I = 0...Q — 1. The Bethe-Ansatz 
may be obtained in the similar way: 

zf = q Q n ^LZJll, / = 1...Q-1. (33) 

m=l,m^l Zl 1 Zm 

The energy is given again by the sum of roots 

E = -q{q-q- x )Y i z h (34) 
i=i 

Inspite of the difference Eqs.( p3| , p4|) must be equivalent to the Eqs.(D,||). 
2. Quadratic Form of Quantum Group Generators 

A limited number of other interesting solvable discrete equations may be obtained from 
a general quadratic form of quantum group generators. Their" classical" version (q — > 1) 
would be differential equations generated by quadratic forms of sli -generators (pO^) These 
differential equations are known in the literature as so-called " quasi exactly soluable " prob- 
lems of quantum mechanics PD| , P2]| . A quadratic form can be considerd as trace of a 
monodromy matrix of an integrable model with nonperiodic boundary conditions 

r = trK + {u)L(su)K^(u)a 2 L\su- 1 )a 2 (35) 

where K± are c-number matrices - solutions of "reflective " Yang-Baxter equations (RYB) 



23| , [24]]. These matrices describe all boundary conditions consistent with integrability 



There is a 3-parametric family of boundary matrices K± which generates a general quadratic 
form of A, B, C, D. Becides, there is a parameter s in (|35D which one can introduce in the 
L-operator (|T7|) preserving integrability. For a particular 

choice of boundary i^-matrices and parameter s r is proportional to the hamiltonian 
(|31~|). Therefore, one can say, that the Hofstadter problem is equivalent to an integrable 



magnet of spin (Q — l)/2 on one site with a proper boundary condition. 
3. q- Analog of Orthogonal Polynomials 

There is an intriguing connection between a wave function fl25]) and q-generalization of 
orthogonal polynomials ||25|| . They satisfy the difference equation (q-analog of differential 
hypergeometrical equation) 
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A{z)P n {q 2 z) + A{z- l )P n {q- 2 z) - {A{z) + A{ Z - x ))P n {z) = {q- 2n - 1)(1 - abcdq 2n ~ 2 )P n {z) (36) 

where = (1 — az)(l — — cz)(l — dz)/((l — z 2 ){l — q 2 z 2 )) and a, b, c, d are parameters 
and n is a degree of the polynomial. Choosing c = — d = q,a = — b = we arrive at the 
equation for the q-Hermite polynomials 

H^\q 2 z) - z 2 H^\q- 2 z) = q- 2n (l - z 2 )H^(z) (37) 

(these are polynomials in z + z~ x of degree n). For an odd Q at n = (Q — l)/2 this yields a 
zero energy solution to eq.([32|) ^ (i?=0) (^) = 2; (0_1)/2 iy ( ( Q ) _ 1 w 2 (2;). 

Another choice c = — d = q, a = —b = q and then the replacement q by q 1 ^ 2 gives the 
q-Legendre equation 

' !/ ^-PW(gz) + iZl^V 1 *) = (g"» + q n+1 )P^\z) (38) 



1-z 2 n w 7 l-z 

Then, comparing with the Eq.(^) we conclude that the zero mode solution is given by the 
q-Legendre polynomial ^^^(iz) = z^^^P^^y^z). 

The almost immidiate and the most interesting task now is to solve the Bethe-Ansatz equa- 
tions in the limit P, Q — ► oo when the flux $/27r is irrational. In all previous examples of 
integrable systems it was always possible to derive an integral equation for a distribution 
function of roots z\. We hope that it would be also possible for this problem and after all 
allows to obtain fractal properties of the spectrum analytically. 
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